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ABSTRACT  OF  OB^CTIVES  AND  ACCOMPLISHMENTS 

A  principal  focus  of  this  research  program  is  the  electronic  structure  aspects  of  nonadiabatic 
processes  related  to  the  stability  and  formation  of  high  energy  density  materials(HEDMs).  Also 
considered  are  nonadiabatic  and  relativistic  effects  that  impact  the  detection  of  molecular  species 
involved  in  HEDM  or  chemcial  laser  development.  These  studies  are  enabled  by  a  unique  system 
of  electronic  structure  algorithms,  the  BROOKLYN  programs,  that  provide  advanced  capabilities 
for  the  study  of  the  electronic  structure  aspects  of  spin-forbidden  and  spin-allowed  electronically 
nonadiabatic  processes.  These  methods,  which  we  believe  define  the  state  of  the  art  in  these 
areas,  have  permitted  us  to  make  significant  contributions  to  the  understanding  of  electronically 
nonadiabatic  processes.  Perhaps  our  most  significant  accomplishement  in  this  grant  period  has 
been  the  development  of  algorithms  for  characterizing  surfaces  of  intersection  of  two  states  of 
distinct  spin-multiplcity  or  two  states  of  the  same  symmetry.  They  complement  our  existing 
capabilities  for  determining  the  relevant  intersurface  couplings,  the  derivative  couplings  that  couple 
states  of  the  same  spin-multiplicity  and  the  spin-orbit  and  dipolar  spin-spin  interactions  that  couple 
states  of  distinct  spin-multiplicity  and  will  be  of  considerable  value  in  studying  the  radiationless 
decay  of  energetic  species.  Of  particular  relevance  in  this  regard  are  our  studies  of  the  stability  of 
the  potential  energetic  species  cyclic  N3F  and  isomers  of  (NO)2.  In  addition  to  their  use  in 
characterizing  the  limits  of  metastability  of  energetic  materials  these  algorithms  will  have  important 
impact  on  our  ability  to  understand  photochemical  processes  involved  in  atmospheric  chemistry 
and  consider  one  of  the  newly  emerging  'hot'  topics  in  molecular  physics  the  importance  of 
Berry's  geometric  phase  in  molecular  processes. 
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SUMMARY 


We  have  made  considerable  progress  in  our  program  to  develop  new  computer  algorithms 
for  the  study  of  spin-conserving  and  spin-nonconserving  electronically  nonadiabatic  processes-. 

We  have  completed  development  of  algorithms  for  characterizing  surfaces  of  intersection  of  two 
states  of  distinct  spin-multiplcity  or  two  states  of  the  same  symmetry.  These  algorithms 
determine  the  surface  of  intersection  of  two  electronic  potential  energy  surfaces  for  which  a  subset 
of  the  intemuclear  parameters  are  kept  fixed  and  the  remainder  are  optimized  to  reduce  the  energy 
of  the  point  of  intersection.  These  algorithms  include  an  automated  procedure  for  constructing  the 
totally  symmetric  portion  of  the  hessian  matrix  using  divided  differences  of  the  appropriate 
analytically  determined  gradients.  They  complement  our  existing  capabilities  for  determining  the 
relevant  intersurface  couplings,  the  derivative  couplings  that  couple  states  of  the  same  spin- 
multiplicity,  and  the  spin-orbit  and  dipolar  spin-spin  interactions  that  couple  states  of  distinct  spin- 
multiplicity,  and  will  be  of  considerable  value  in  studying  the  radiationless  decay  of  energetic 
species.  In  addition  to  their  use  in  characterizing  the  limits  of  metastability  of  energetic  materials 
these  algorithms  will  have  important  impact  on  our  ability  to  understand  photochemical  processes 
involved  in  atmospheric  chemistry  and  consider  one  of  the  newly  emerging  'hot'  topics  in 
molecular  physics  the  importance  of  Berry's  geometric  phase  in  molecular  processes. 

Our  need  to  estimate  the  lifetime  of  metastable  species  has  lead  us  to  develop  a 
nonperturbative  computational  approach  to  treat  resonances  attributable  to  nonadiabatic 
radiationless  decay  into  several  electronically  distinct  channels.  Expressions  for  both  the  radiative 
and  radiationless  decay  rates  were  developed  and  an  alternative  derivation  of  Fano's  classic 
multichannel  line  shape  formula  was  obtained. 

Previously  we  investigated  the  stability  of  the  potential  energetic  materials,  a-N202.  We 
showed  that  a-N202,  stable  with  respect  to  spin-allowed  decay,  is  in  fact  rapidly  predissociated 
to  N2O  +  0(3p)  via  a  spin-forbidden  channel.  We  have  continued  our  studies  along  these  lines 
completing,  in  collaboration  with  M.  Gordon  at  Iowa  State  University,  a  study  of  the  spin- 
forbidden  predissociation  of  the  potential  HEDM  cyclic  N3F,  NsFCUA')  N2(X'E^  +  NFCX^E- 
).  It  was  concluded  that  this  system  is  stable  with  respect  to  spin-forbidden  decay  to  the  N2  + 
NF(X3e-)  asymptote 


The  spectroscopy  of  small  boron  containing  compounds  is  of  consequence  owing  to  the 
importance  of  boron  as  a  fuel  additive.  Brazier  of  the  Phillips  Laboratory,  has  pointed  out 
existence  of  significant  nonadiabatic  effects  in  the  l,23ng  states  of  B2  that  complicate  the  analysis 
of  the  B2  spectra,  particularly  the  vibrational  level  dependence  of  the  fine  structure  interaction. 
Although  the  electronic  states  of  B2  had  been  the  object  of  several  high  quality  theoretical  studies 
the  nonadiabatic  interactions  had  not  been  characterized.  As  part  of  this  grant  the  vibrational  level 
dependence  of  the  fine  structure  splitting  was  determined. 

The  metastable  a^ff  state  of  BH  is  involved  in  a  potentital  energy  transfer  chemical  laser, 
suggested  by  D.  Benard,  based  on  the  BH(Aln)  ->  BH(XiZ+)  electronic  transition.  The  radiative 
lifetime  of  this  species,  attributable  to  the  spin-forbidden  transition  a^fl  —>  is  an  important 
factor  in  the  design  of  such  a  laser  system.  The  lifetimes  of  the  individual  fine  structure  levels 
(a^n,  V,  N,  Fj,  e/f)  were  determined  in  a  pure  Hund's  case  (b)  approximation.  In  all  cases  the 
(a^n,  V,  N,  Fi,  e/f)  levels  are  long-lived  with  lifetimes  longer  than  10  s  for  N  =  0  and  v  =  0,1, 2,3. 

Because  of  the  potential  role  of  BH/a^ri)  in  the  above  noted  chemical  laser  system 
quantitative  detection  schemes  are  desirable.  Laboratory  detection  of  the  a^fl  state  is  based  on 
b^Z~  — >  a^n  laser  induced  fluorescence.  In  order  for  this  transition  to  be  used  to  characterize  the 
a^n  state  in  a  quantitative  manner  the  radiative  and  radiationless  decay  rates  of  BH/b^Z")  are 
required.  To  this  end  the  rates  for  the  spin-allowed  radiative  transitions  (b^Z",  v,  J)  ^  a^II  and 
the  spin-orbit  induced  radiationless  decay  process  b^Z-  ^  l3i;+  were  determined.  It  was  found 
that  predissociation  of  BH(b^Z~,  v)  is  negligible  up  to  v=2,  where  measureable  predissociation 
exists.  This  research  was  performed  in  collaboration  with  P.  J  .Dagdigian's  group  at  Johns 
Hopkins.  This  experimental-theoretical  collaboration  was  served  to  enhance  considerably  the 
impact  of  this  study. 

TECHNICAL  REPORT 

The  thrust  of  our  research  program  is  to  develop  computational  techniques  to  determine  the 
stability  with  respect  to  radiationless  decay,  of  potential  high  energy  density  materials.  These 
techniques  are  then  to  be  applied  to  the  study  of  chemically  relevant  systems.  We  are  also 
concerned  with  spectroscopic  analyses  relevant  to  the  characterzation  and  detection  of  energetic 
species.  During  the  present  performance  period  we  have  produced  new  tools  for  the  study  of 
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nonadiabatic  processes.  In  addition  to  their  use  in  studying  to  the  stability  of  energetic  materials 
these  algorithms  have  potentially  broad  impact,  including  applications  in  atmospheric 
photochemistry  and  as  a  tool  to  consider  one  of  the  newly  emerging  'hot'  topics  in  molecular 
physics  the  importance  of  Berry's  geometric  phase  in  molecular  processes.  Below  we  discuss  the 
most  significant  components  of  our  current  research  effort  including  (lA)  development  of  a 
computational  methodology  to  characterize  surfaces  of  intersection  of  two  states  of  either,  distinct 
spin-multiplicity  or,  the  same  overall  symmetry,  (IIA)  the  use  of  this  methodology  to  consider,  in 
collaboration  with  M.  Gordon  at  Iowa  State  University,  the  possibility  of  spin-forbidden 
radiationless  decay  in  the  potential  energetic  material  cyclic  N3F,  (IIC)  the  radiative  lifetime  of 
BHCa^II)  an  intermediate  in  a  chemical  laser  system  proposed  by  D.  Benard,  and  (IID)  the  radiative 
and  radiationless  decay  rates  of  BHCa^Z-)  required  for  the  quantitative  detection  of  BHCa^Tl)  using 
the  b^E”  — >  a^n  transition  in  laser  induced  fluorescence. 


I.  New  Tools  for  Characterizing  Nonadiabatic  Processes 

A.  Characterization  of  Surfaces  of  Intersection  of  Two  States  of  Distinct  Multiplicity  or  Two 
States  of  the  Same  Symmetry  Using  Analytic  Derivative  Techniques 

When  considering  the  stability  of  an  energetic  species  corresponding  to  a  local  minimum  on 

the  ground  potential  energy  surface  one  must  consider  the  possibility  of  radiationless  decay 

induced  by  intersections  with  nominally  excited  potential  energy  surfaces.  The  straightforward 

procedure  that  has  been  used  in  the  past  to  characterize  surfaces  of  intersection  initially  determines 

the  potential  energy  surfaces  in  question  and  then  characterizes  their  intersections.  We  have 

developed ’’2  an  algorithm  that  directly  determines,  that  is  without  prior  determination  of  the 

potential  energy  surfaces  in  question,  surfaces  of  intersection  of  two  states  of  distinct  multiplicity 

or  two  states  of  the  same  symmetry  using  analytic  derivative  techniques. 


The  algorithm  determines  the  solution  of  the  following  Newton-Raphson  equations: 

■q'^(R,U)  g"(R)  h^^(R)  k(R)ir5Rl  rg'(R)  +  ^,g^^(R)  +  ^X(R)  +  ^;"^_A,k'(R)l  1-la 
g^-'(R)^  0  0  0  AE^CR)  1-lb 

h'^(R)'  0  0  0  <5^2  ~  0  1.1c 

k(R)'  0^  0^  0  l.ld 

where  ^  and  X  are  Lagrange  multipliers,  AE,°  (R)  =  E^{R)  -  eJ(R),  the  K(R)  are  a  set  of 
geometrical  constraints,  5R  =  R'-  R,  6X=X'  -  X,  6^  =  ^'  -  4,  g^(R)  =  , 
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Sa  (R)  =  — - >  "a  (R)  =  ^  ^  >  ^a(R)  =  — :3 —  and  Q^-’  is  a  matrix  of  second 

derivative  terms.  Implicit  in  the  design  of  the  algorithm  based  on  this  system  of  equations  is  the 
assumption  that  the  gradients /z"(R),  g"(R)  and  g'(R)  are  readily  available.  This  is 
accomplished  through  the  use  of  analytic  gradient  techniques. 

Eq.  1.1a  achieves  the  energy  minimization.  Eqs  l.lb,c  enforce  the  requirement  Ei(R)  = 
Ej(R)  while  eq.  1.  Id  enforces  the  geometrical  constraints,  K(R)  =  0.  Efficient  energy 
minimization  requires  the  availability  of  the  Hessian  Q^J(R,  A,).  The  R,  ^  contribution,  QIJ(R, 
0)  is  currently  obtained  by  divided  difference 


QiJ(R,  t  0)  =  [(g'  +  +4h")(R  +  er)  -  (g^  +  +^3h'^)(R)]/e  1.2 

where  1“  is  a  unit  vector  in  the  direction  R“.  The  X  contribution  to  Q^J(R,  A,)  is  readily 
evaluated  analytically.  The  1“  are  usually  taken  to  be  atom  centered  cartesian  coordinates,  Wy, 
where  W  =  X,Y,Z,  and  y  =  1-  M  denotes  the  M  atoms.  However  when  spatial  symmetry 
constraints  are  imposed  use  of  eq.  1.2  for  1“  =  Wyis  computationally  inconvenient.  These 
displacements  do  not  preserve  the  point  group  symmetry',  so  that  symmetry  cannot  be  exploited  in 
the  evaluation  of  the  electronic  wavefunction.  In  addition  when  spatial  symmetry  is  exploited  the 
number  of  degrees  of  freedom  that  must  be  considered  in  the  evaluation  of  eq.  1 .2,  a  time 
consuming  step,  is  reduced.  An  automated  procedure  has  been  implemented  to  evaluate  eq.  1.2 
with  respect  to  the  totally  symmetric  internal  coordinates.  This  algorithm  has  been  used  to 
considerable  advantage  in  the  studies  of  cyclic  N3F  discussed  below. 

B .  Nonadiabatic  Radiationless  Decay  of  Resonances 

To  faciliate  our  studies  of  radiationless  decay  diatomic  systems  a  nonperturbative 
computational  approach  was  developed  to  treat  resonances  attributable  to  nonadiabatic  radiationless 
decay  into  several  electronically  distinct  channels.  Expressions  for  both  the  radiative  and 
radiationless  decay  rates  were  developed  and  an  alternative  derivation  of  Fano's  classic 
multichannel  line  shape  formula  was  obtained.  Although  the  approach  is  nonperturbative  the 
linewidth  is  expressed  in  a  Golden  mle  type  formula.  The  resonance  energy  is  obtained  from  the 
iterative  solution  of  an  eigenvalue  problem  in  the  bound  state  space  avoiding  the  tedious 
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characterization  of  the  5-matrix  as  a  function  of  scattering  energy.  These  attributes  enable  efficient 
determination  both  narrow  and  broad  linewidths  and  large  resonance  energy  shifts. 

II.  Applications 

A.  Spin-forbidden  Radiationless  Decay  of  Energetic  Materials 

To  address  the  question  of  spin-forbidden  radiationless  decay  of  an  energetic  material  in  a 
singlet  electronic  state  the  topology  of  the  singlet  -  triplet  surface  of  intersection  is  required.  The 
range  of  situations  that  can  arise  is  illustrated  in  the  following  figure; 


where  intersection  i  lies  inside  the  barrier  to  spin-allowed  dissociaton,  intersection  3  lies  behind 
that  barrier  and  intersection  2  represents  an  intermediate  case.  In  situation  3(1)  the  species  will 
(not)  be  sufficiently  long  lived  to  be  a  potential  HEDM  species.  In  situation  2  careful  consideration 
of  the  nonadiabatic  dynamics  will  be  required. 

1.  N3F(1A’)  N2(Xli:p  +  NF(X3Z-) 

As  part  of  an  ongoing  collaboration  with  Gordon's  group  at  Iowa  State  University  the  ^A' 
and  3A"  potential  energy  surfaces  for  the  N3F  molecule  were  studied  using  multi-configurational 
wave  functions.^  Two  new  isomers  were  found,  one  on  the  singlet  (lA')  and  one  on  the  triplet 
(3A")  surface.  Both  isomers  have  a  3-membered  cyclic  structure  and  Cg  symmetry.  The  singlet 
isomer  is  pictured  below 


6 


The  singlet  cyclic  isomer  is  endoergic  relative  to  the  open  fluorine  azide  by  15-17  kcal/mol.  Its 
kinetic  stability  is  close  to  the  stability  of  the  open  isomer:  the  barrier  separating  the  cyclic  isomer 
from  the  dissociation  products  N2(XlIp+NF  (al  A)  is  about  13-17  kcal/mol  and  is  lower  than  the 
barrier  to  isomerization.  The  triplet  cyclic  isomer  is  much  higher  in  energy  (about  70  kcal/mol), 
with  a  barrier  to  dissociation  to  N2(X1S^-i-NF(X3e-)  on  the  order  of  15  kcal/mol  but  can  decay 
radiatively  to  the  lower  singlet  surface.  Crossings  of  the  ^A'  and  the  ^A"  surfaces  may  allow  the 
cyclic  singlet  isomer  to  predissociate  to  the  ground  state  products,  N2(XlZ^4-NF(X3z-).  It  was 
shown,  however,  that  the  singlet-triplet  surface  of  intersection  lies  'behind'  (intersection  3)  the 
barrier  to  singlet  decomposition,  so  that  spin-forbidden  predissociation  will  not  preclude  detection 
of  cyclic  N3F. 

(ii).  N202(lAig)(D2h)^  NO(X2n)  +  NOfX^n) 

The  geometric  isomers  of  (NO)2  have  long  been  of  interest  to  the  HEDM  program  since 
Michels  and  Montgomery^  first  suggested  a-N202  (NNOO)  as  a  possible  energetic  material. 
Subsequently  we  showed  that  this  isomer  is  unstable  with  respect  to  spin-forbidden  radationless 
decay  to  N2O  4-  0(3p)  (intersection  1 )?  Currently  we,  again  in  collaboration  with  Chaban  and 
Gordon,  are  considering  the  stability  of  the  D2h  form  of  N2O2,  (^Aig).  This  species  is  metastable 
with  respect  to  two  NO(X2n)  and  has  a  barrier  to  spin-allowed  decay  of  approximately  (MCSCF 
level)  47  kcal/mol.  In  this  instance  preliminary  calculations  suggest  that  a  singlet-  triplet  crossing 
occurs  inside  the  barrier  to  spin-allowed  predissociation.  Calculations  are  currently  in  progress  to 
characterize  the  singlet-triplet  surface  of  intersection. 

B.  Nonadiabatic  Effects  in  the  Spectra  of  HEDM  Species 

(i)  The  Spectra  of  the  l,23ng  States  of  B2 

Our  interest  in  B2  has  been  stimulated  the  potential  of  boron  containing  compounds,  as 
high  energy  density  propellants.^  Recently  observation  of  five  band  systems  of  B2^  has  been 
reported,  including  two  bands  involving  the  l,23ng  states,  the  23nu  -  l3ng  and  the  23ng  - 
A^riu  bands.  These  studies  benefitted  from  several  high  quality  theoretical  studies  of  this 
homonuclear  diatomic  molecule.^‘*2  Despite  the  availability  of  high  quality  ab  initio  potential 
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energy  curves  and  interstate  electronic  transition  moments,  detailed  spectral  analyses  are  hampered 
by  the  presence  of  nonadiabatic  perturbations  in  the  spectra.  Of  paiticular  concern  to  our  research 
are  the  l,23ng  states.  Brazier  and  Garrick^  observed  significant  deviations  from  the  expected 
isotopomer  relations  in  the  (23ng,v  =  0)  -  (A3nu,v=  0)  band,  deviations  they  attributed  to 
perturbations  in  the  23ng  state.  On  the  basis  of  previous  theoretical  work*^  these  perturbations  are 
likely  to  result  from  an  avoided  crossing  with  the  Pflg  state.  The  avoided  crossing  of  the  l,23ng 
states  is  expected  to  significantly  effect  the  vibrational  level  dependence  of  the  fine  structure 
splitting  constant  in  both  these  states.  This  quantity  is  of  practical  importance  as  the  inference  of 
the  fine  structure  splitting  constant  was  key  to  a  recent  assignment^  of  the  (23nu,v=  0)  -  (l3ng,v 
=  0)  band. 

For  these  reasons  nonadiabatic  effects  and  fine  structure  splittings  in  the  l,23ng  states  of 
B2  were  studied.  Adiabatic  electronic  wavefunctions  were  computed  at  the 

multireference  configuration  interaction  level.  The  interstate  derivative  couplings /-^(R)  = 

,  7=1,  7=  2,  were  computed  and  used  to  construct  'rigorous' 
diabatic  electronic  states.  The  spin-orbit  interactions  responsible  for  the  fine  structure  splitting  in 
the  l,23ng  states  were  determined.  These  electronic  structure  data  were  used  to  determine  the 
vibrational  structure  of  the  l,23ng  states  in  both  coupled  adiabatic  and  coupled  diabatic  states 
bases.  As  a  result  of  an  avoided  crossing  of  the  l,23ng  states  the  fine  structure  splitting  constants 
exhibit  a  marked  geometry  dependence.  The  fine  structure  splitting  constants  for  the  (l^flg,  v  =  0) 
and  (23ng ,  v  =  0)  levels  were  found  to  be  -4.2(-4.4)^,  1.37(1.54)'3  cm-i,  respectively  in  good 
agreement  with  a  recent  experimental  determinations  given  parenthetically. 

(ii)  Predissociation  of  the  l,23ng  electronic  states  of  B2  -  Work  In  Progress 
Brazier  and  Garrick' 3  observed  a  sharp  cut  off  in  the  vibrational  spectrum  of  the  23ng  state 
of  B2  for  V  >  1 .  One  of  the  goals  of  our  current  research  effort  is  to  extend  the  above  study  to 
provide  a  theoretical  underpining  for  this  observation. 

To  address  this  question  the  resonance  theory  of  subsection  IB  is  being  used  in  conjunction 
with  our  capabilities  to  determine  the  derivative  couplings  /'''(R)  noted  above  .  In 
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contradistinction  to  the  above  study  however,  here  adiabatic  electronic  wavefunctions  for  three 
3ng  states  are  required.  This  has  been  accomplished  at  the  multireference  configuration  interaction 
level.  The  associated/-^(R)  were  computed  and  used  to  construct  'rigorous'  diabatic  electronic 
states.  These  electronic  structure  data  are  being  used  to  understand  the  mechanism  of  the 
predissociation  of  the  {2^Ylg,  v  =  2)  level  which  in  the  adiabatic  state  approximation  is  a  purely 
bound  state  but  is  not  observed  spectroscopically.'^ 

C.  Radiative  Decay  of  Electronically  Excited  Metastable  Species 

(i)  The  Lifetime  of  BH(a3n) 

Benard  et.  have  reported  initial  efforts  toward  the  development  of  a  chemical  laser 

system  based  on  the  A'll  —  X'S"*"  transition  in  BH.  BH(A'n)  is  produced  by  energy  transfer  from 
the  metastable  a'A  state  of  NF,  Te  =  1 1442  cm-l,'^  as  follows,  BH(X'2:+)  +  NFfa^A)  -» 
BHfa^ri)  +  NF(X3z~)  and  BHfa^ll)  +  NF(a^A)  — >  BH(A'n)  +  NF(X3Z“).  This  work  has 
stimulated  our  interest  in  the  stability  and  detection(see  below)  of  BHCa^n).  To  this  end  the  spin- 
forbidden  radiative  decay  process  BHCa^fl)  BH(X'I+)  was  studied'’^  using  extended 
correlation  consistent  basis  sets,  second  order  configuration  interaction  wavefunctions  and  the 
Breit-Pauli  spin-orbit  approximation.'*  The  two  unique  electronic  transition  moments,  |ix  and  |i||, 
connecting  the  a^fl  and  X'X+  nonrotating  molecule  states  were  determined.  The  lifetimes  of  the 
individual  fine  structure  levels  (a^n,  v,  N,  Fi,  e/f)  were  determined  in  a  pure  Fund's  case  (b) 
approximation.  In  all  cases  the  (a^FI,  v,  N,  Fj,  e/f)  levels  are  long-lived  with  lifetimes  longer  than 
10  s  for  N  =  0  and  v  =  0,1,2,3. 

D.  Spectroscopic  Identification  of  Chemical  Species 

(i)  Radiationless  decay  of  BH/b^Z") 

Laboratory  detection  of  BH(a3n)  is  based  on  b^Z-  ^  a^II  laser  induced  fluorescence.'^ 

In  order  for  this  transition  to  be  used  to  characterize  the  a^n  state  in  a  quantitative  manner  the 
radiative  and  radiationless  decay  rates  of  BHCb^Z")  are  required.  To  this  end  rates  for  the  spin- 
allowed  radiative  transitions  (b^Z-,  v,  J)  a^n  (Ref  17  )  and  the  spin-orbit  induced  radiationless 
decay  process  b^Z"  (Ref.  20 )  were  determined.  The  spin-orbit  coupling  was  determined 
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within  the  Breit-Pauli  approximation  using  MCSCF/CI  wavefunctions.  The  radiationless  decay 
lifetimes  of  the  (b^Z",v,N,  Fi)  fine  structure  levels  were  determined  in  a  Hund's  case  (b) 
approximation.  The  (b3Z‘,v  =  0-2,  N  =  0-15)  levels  are  not  significantly  predissociated. 
However  a  spin-orbit  induced  perturbation  b^Z'  ~  l3x+  of  magnitude  approximately  0.2  cm’l 
within  the  crossing  region  of  the  h^Ir  and  l3l+  states  results  in  appreciable  predissociation  of 
the  v'  =  3  (v'  =  2)  level  in  BH  for  low  N  (high  N). 

These  theoretical  predictions  were  compared  with  lifetimes  of  rotational/fine-structure  levels 
of  electronically  excited  BH(Z53l-),  in  vibrational  levels  v'  =  0  -  4,  determined  from  fluorescence 
decay  waveforms  with  laser  excitation  on  isolated  rotational  lines  in  the  b^IT  -  a^U  Av  =  0 
sequence.^’  The  experimental  and  theoretical  lifetimes  agree  very  well.  This  comparison  confirms 
that  the  low  rotational  levels  for  v'  =  0  -  2  decay  purely  radiatively,  with  a  rate  decreasing  as 
v'  increases.  These  rates  are  virtually  independent  of  the  fine-structure  level.  However  the 
lifetimes  of  the  higher  vibrational  levels  exhibit  significant  fine-structure  level  dependence,  with  the 
shortest  values  obtained  for  the  F2  levels.  This  behavior  was  shown  to  be  attributable  to 
predissociation  induced  by  b^lT  ~  PZ'*'  spin-orbit  coupling. 

Ill  Future  Directions 

(i)  Prediction  of  MCD  Spectra 

The  characterization  of  the  local  stmcture  of  first  row  elements  in  hydrogen  matrices  is  an 
important  problem  in  the  study  of  cryogenic  propellants.  Magnetic  circular  dichroism  (MCD)  is 
one  of  the  tools  currently  being  used  to  investigate  the  problem.  We  are  in  the  process  of 
developing  algorithms  to  calculate  MCD  spectra  from  first  principles.  As  discussed  in  detail  on  our 
renewal  proposal,  our  principle  methodological  advance  involves  the  key  B-term.  We  will  replace 
the  traditional,  but  computationally  intractable,  sum  over  states  ( the  K-sum) 

X  ,X,K^J 

I™  Z(n.  )  X  (n  N'*'*)  ■  ('*'*  )> 

K^A,X,X' 

with  a  direct  solution  of  a  system  of  linear  equations 


10 


=  imX(>F'.  H'f'OxK  H'^'T} 

A.A’  A.A' 


where 


This  work  will  support  an  AFOSR  funded  experimental  program  in  Professor  Kenney's 
laboratory. 

(ii)  Pathways  to  energetic  species 

Conical  intersections  of  two  states  of  the  same  symmetry  provide  efficient  pathways  for 
electronically  excited  states  to  decay  radiationlessly  to  the  ground  electronic  state.  The  barriers  that 
'protect'  metstable  structures  suggest  the  nearbye  presence  of  conical  intersections  with  excited 
states.  These  conical  intersections,  funnels,  may  provide  pathways  to  the  metastable  structures 
(stabilized  by  collisions)  that  represent  HEDM  species.  While  certainly  not  a  viable  synthetic  route 
for  large  quantities  of  such  a  material  these  pathways  would  allow  for  the  experimental  study  of 
such  systems.  The  algortihms  describe  above  will  be  used  to  consider  this  possibility. 
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February  10,  1995 

Nonadiabatic  Chemistry:  New  Computational  Challenges 
presented  at:  Iowa  State  University 
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Department  of  Chemistry, 
March  31,  1995 


Nonadiabatic  Chemistry:  NewTheoretical  Challenges 

presented  at:  Cal  Tech 

Department  of  Chemistry, 

May  2,  1995 

Nonadiabatic  Chemistry:  NewTheoretical  Challenges 

presented  at:  State  University  of  Leiden 
Gorlaeus  Laboratory 
Leiden,  The  Netherlands 
July  6,  1995 

TECHNOLOGY  TRANSFER 

Radiative  and  Radiationless  decay  rates  of  state  of  BH 

D.  J.  Benard,  Ellen  Boehmer 

Science  Center,  Rockwell  International  Corp. 

and 

C.  R.  Brazier  and  P.  G.  Carrick 
Phillips  Laboratory,  Edwards  AFB 

transfer:  To  quantify  conditions  in  a  laser  system  LIE  may  be  used  a  probe.  In  order  to  be  a 
quantitative  probe  the  radiative  AND  radiationless  decay  rates  of  the  laser  excited  state  must  be 
known.  Lisa  Pederson  calculated  the  radiative  decay  rate  for  the  b^Z”  ^  a^II  transition  (Pub.  85) 
and  the  decay  rate  for  the  radiationless  transition  h^Ir  ~  \^1+  (Pub.  88)  in  BH.  These  results 
were  communicated  by  Ms  Pederson  to  E.  Boehmer.  This  information  is  needed  to  determine 
BH(a3n,  v)  level  populations  by  LIE.  BHla^H)  is  an  energetic  intermediate  in  a  potential  laser 
system  being  considered  by  Rockwell. 

We  have  communicated  the  results  of  our  combined  experimental/theoretical  treatment  [Pubs.  88 
and  94]  of  the  decay  of  the  b^Z-  state  to  C.  Brazier  of  the  Phillips  Laboratory,  Edwards  APB. 
This  will  interaction  will  be  useful  in  understanding  Brazier's  measured  spectra  of  BH  b^Z- 

Comments:  The  capability  to  determine  spin-allowed  radiative  decay  rates  is  routine  in  groups 
specializing  in  electronic  structure  theory.  However  the  capability  to  quantify  the  spin-orbit 
induced  radiationless  decay  process  b^Z-  ~  l3z+  is  not.  Our  capabities  define  the  state  of  the 
state  in  this  area  and  have  important  applications  in  various  processes  of  interest  to  the  Air  Porce  . 

INTERACTIONS 


C.  R.  Brazier  and  P.  G.  Carrick 
Phillips  Laboratory,  Edwards  APB 
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We  have  performed  a  theoretical  analysis  of  the  vibrational  level  dependence  of  the  fine 
structure  splitting  constant  in  the  l,23ng  states  of  B2.  This  analysis  required  the  determination  of 
nonadiabatic  interactions,  derivative  couplings,  between  the  l,23ng  adiabatic  states  as  well  is  the 
spin-orbit  interactions.  We  are  currently  extending  this  work  to  include  predissociation  of  the 
l,23ng  states.  This  work  [Pub.  83  and  work  in  progress]  was/will  be  communcated  to  Brazier 
and  Garrick  and  supports  their  spectroscopic  studies  of  these  states. 

Comments:  The  capability  to  determine  the  requisite  potential  energy  curves  is  routine  in  groups 
specializing  in  electronic  structure  theory.  However  the  capability  to  quantify  the  nonadiabatic 
(derivative)  couplings  is  not.  Again  our  capabilities  define  the  state  of  the  art  in  this  area. 
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